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Abstract 

y— I ■ 

, The Ehlers-Ellis 1+3 formulation of covariant hydrodynamics, when sup- 

qq | plemented with covariant radiative transport theory, gives an exact, physically 

CN ■ transparent description of the physics of the cosmic microwave background ra- 

diation (CMB). Linearisation around a Friedmann-Robertson- Walker (FRW) 
universe provides a very direct and seamless route through to the linear, gauge- 
invariant perturbation equations for scalar, vector and tensor modes in an 
almost-FRW model. In this contribution we review covariant radiative trans- 
port theory and its application to the perturbative method for calculating and 
understanding the anisotropy of the CMB. Particular emphasis is placed on 
the inclusion of polarization in a fully covariant manner. With this inclusion, 
the covariant perturbative approach offers a complete description of linearised 
CMB physics in an almost-FRW universe. 

K 

1 Introduction 

Given the potential impact of observations of the cosmic microwave background 
radiation on our understanding of the large scale properties of the universe, there 
is a strong case for developing a physically transparent and portable formalism for 
the study of CMB physics. Ideally, such a framework should allow one to give: (i) 
an exact, but physically illuminating description of all the relevant physics; (ii) a 
simple, gauge-invariant linear perturbation theory around a variety of background 
models. Along with many other workers in this field 0-0, we have found the 1+3 
covariant formulation of hydrodynamics 0, |J and radiative transport theory |T0| , 11 
to be well suited to this task 
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In this contribution we begin by reviewing the exact 1+3 covariant theory of 
polarized radiative transport in the presence of Thomson scattering (the dominant 
scattering mechanism for CMB photons over the epoch of interest). This work builds 
on the covariant approach to kinetic theory developed in |10], [IH| . The angular de- 
pendence of the intensity and polarization are conveniently described in terms of 
projected symmetric trace free (PSTF) tensors. The multipole equations govern- 
ing the dynamics of these radiation variables are given in linearised form around 
an almost-FRW universe, and generic features of their solution are discussed qual- 
itatively. By further decomposing the linearised multipole equations in harmonic 
functions, we arrive at a complete set of gauge-invariant perturbation equations in 
a form which is convenient for numerical implementation, but which still preserves 
the physical transparency of the 1+3 covariant approach. We present the mode- 
expanded perturbation equations for scalar (density clumping) modes, and discuss 
briefly their numerical implementation and some numerical results. 

We employ standard general relativity and adopt a (H ) signature for 

the metric g a b. Our conventions for the Riemann and Ricci tensors are fixed by 
[V fl , Vfc]w c = R a bc d Ud and R ac = R a bc- The spacetime alternating tensor is de- 
noted r] a bcd- Round brackets denote symmetrisation on the enclosed indices, square 
brackets denote antisymmetrisation, and angle brackets denote the PSTF part (the 
projection is with h a b = g a b — u a Ub, where u a is the fundamental velocity of the 1+3 
approach). Overdots denote the action of w a V a , and the totally projected covariant 
"derivative" is denoted by D a . 



2 Covariant Radiative Transport Theory 

In the 1+3 covariant approach to cosmology (see fl2| for a recent review), the cos- 



mological model is described in terms of geometric variables that have a direct in- 
terpretation in terms of the observations of a preferred set of observers who are 
comoving with the fundamental velocity field u a . We shall refer to the choice of u a 
as a frame choice^. We consider a diffuse radiation field (the CMB) observed from 
the u a frame. The Stokes parameters of the radiation in some solid angle dVL and 
some photon energy range dE, along the unit projected vector e a , can be combined 
into a second-rank polarization tensor P a b(E, e) which is projected relative to u a and 
e a on both indices. The non-vanishing components of P a b(E,e) on an orthonormal 
tetrad with the O-direction aligned with u a and the 3-direction with e a are defined 

lr The issue of how to specify u a in an appropriate physical manner is discussed in [T^ |. 
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in terms of the Stokes parameters via: 



P n = ±[I(E,e) + Q(E,e)) 
P 12 = l[U(E,e) + V(E,e)] 
P 21 = ±[U{E,e)-V{E,e)\ 

P 22 = §[!(£, e)-Q(£,e)]. (1) 



Our notation for the Stokes parameters follows ||13f| . The transformation properties 
of the Stokes parameters under rotation of the 1- and 2-directions ensures that 
P ab (E,e), as defined, is tensor-valued. Introducing a screen projection tensor H ab , 

T-tab = Kb + e a e b , (2) 

the polarization tensor can be written in irreducible form 

P ab (E, e) = -\I(E, e)H ab + P ab (E, e) - \V(E, e)e abc e c , (3) 

where e abc = r] abcd u d and V a b{E, e) is a PSTF tensor which is orthogonal to e a : 

V ab (E,e)=H c a V cb (E,e). (4) 

2.1 The Multipole Decomposition 

For many applications it is convenient to describe the angular dependence of the 
polarization tensor in terms of PSTF tensor-valued multipole moments. Since I(E, e) 
and V(E, e) are scalar fields on the sphere they can be expanded as follows: 

J(£,e) = E J ^) eA ' ( 5 ) 

1=0 

oo 

V(E,e) = Y;V Al (E)e A >, (6) 



1=0 



where Iai(E) and Va^E) are energy-dependent PSTF tensors, and we have intro- 
duced the index notation Ai = ai . . . ai with e Al = e ai . . . e ai . The expansions of 
I(E, e) and V(E, e) can be inverted to determine the multipoles, for example 



I Al (E) = ^Jdni(E,e)e {Alh with A^ ^^J . (7) 



For the linear polarization tensor V ab {E,e) we must expand in tensor spherical 
harmonics that are projected perpendicular to e a (tensors projected perpendicular 
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to e a , such as V a b(E, e), are denoted transverse in WW). It is most convenient to use 



pure-spin PSTF harmonics ||14|| , so that 

00 oo 

P aia2 (E,e) = ^[£: aia2B ^ 2 e^-] TT + ^[e cl e ClC2 {ai S a2)c2B; _ 2 e^-] TT . (8) 

1=2 1=2 

Coordinate-dependent versions of this expansion were first discussed in the context 
of all-sky CMB polarization data in [|l|, |16fl . In equation (|) we have left the energy 
dependence of the PSTF tensors S Al (E) and B Al (E) implicit. The superscript tt 
denotes the transverse traceless part: 

[A ab ] TT = n:m c b 2 A clC2 - \n ab n c ^A clC2 , (9) 

where A ab is a second-rank tensor. The two summations in equation (H) define the 
electric and magnetic parts of V a b{E,e). The l-th term in the electric part has 
parity (—1)', while the l-th term in the magnetic part has parity (— 1) . The 
electric and magnetic multipoles can be determined from the polarization tensor 
with the following inversion formulae: 

£ Al (E) = ^ (l+l)^ 1 ^) J Me^P^&e), (10) 

-B Al {E) = + % + J dne b e bc {ai e Al _ 2 V ai _ l)c (E,e), (11) 

which follow from the orthogonality of the eaj. We shall also make use of the 
energy-integrated multipoles, I Al , £-a v B Al and V Al , where, for example, 

I Al = A t J dEI Al (E). (12) 

It follows that the three lowest energy-integrated multipoles of I(E, e) give the ra- 
diation energy density, momentum density and anisotropic stress respectively: 

1 = P {1 \ Ia = ti\ h b = ^ b ] . (13) 

The scale dependence of intensity and polarization correlations on the CMB sky 
are described by their angular power spectra. In a statistically isotropic ensemble 
of CMB skies, the power spectrum of intensity anisotropies Cf 1 is defined by the 
ensemble average || 

(j)\W B *) = Wf%,h™, (14) 

where h/^l = h^...h b A. In the limit of almost-isotropic radiation, the I Al are 
related to the gauge-invariant bolometric temperature perturbation (from the all- 
sky mean) <5r(e) via J2| 

00 -i 

^^EaV 1 , (15) 
i=i 1 
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from which we recover the usual Legendre expansion of the temperature correlation 
function: 



oo 



(8 T (e)6 T (e')) = ^ i^^Cf P,(- e V a ). (16) 
1=1 

Power spectra for the polarization and the cross-correlation between electric po- 
larization and intensity anisotropies can be defined in a similar manner to Cf 1 . 
However, for consistency with previous approaches based on an expansion in spin- 
weighted harmonics |16|], with the normalisation conventions adopted here it is nec- 



essary to include an additional factor of — l)/(/ + 1)(/ + 2)] in the definition of 
Ci for each factor of the polarization. For example, 

(l + 2)(l + 

Note there cannot be any correlation between Ba 1 and either of S^ t or in a parity 
symmetric ensemble. As noted in |16||, the definitions of the polarization power 



spectra adopted here are independent of the choice of polarization basis vectors, 
which has significant advantages over correlating the Stokes parameters directly. 



2.2 The Boltzmann Equation 

The evolution of the polarization tensor P ab (E, e) along the path of the radiation in 
phase space is described by the (collisional) Boltzmann equation 

C[E- 3 P ab (E,e)] = K ab (E,e). (18) 

The appropriate Liouville operator £ for the transverse tensor P ab (E, e) acts accord- 
ing to 

CP ab (E, e) = Him?p d V d P cl ci(E, e), (19) 

where the derivative is along the photon path {x a (X),p a (X)} in phase space. Here, 
p a = E(u a + e a ) is the photon momentum and A is an affine parameter with p a = 
dx a /d\. The scattering tensor K ab (E,e) on the right-hand side of equation (|18|) is 
transverse. Its exact form for Thomson scattering is given in Section p. 2. 2 . 

The irreducible decomposition of P ab (E,e) into I(E,e), V(E,e) and V ab (E,e), 
given as equation (|3|), is preserved by the action of the Liouville operator, so that 

C[E- 3 P ab (E,e)} = -^[E- 3 I(E,e)]H ab + C[E- 3 V ab (E,e)} 

-\^[E- 3 V{E,e)]e abc e c . (20) 

In the absence of scattering (K ab (E, e) = 0), we recover the well-known results that 
the occupation number I(E,e)/E 3 , degree of circular polarization V(E,e)/I(E,e), 
and degree of linear polarization y/[2V a b(E, e)V ab (E, e)]/I(E, e) are independently 
conserved along the path of the radiation. 
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2.2.1 Transformation Laws Under Changes of Frame 

The 1+3 covariant approach considers observations made from the viewpoint of a 
set of observers moving with the fundamental velocity u a . Often, it is useful to have 
available the transformation laws relating these observations to those made with 
respect to a different velocity field u a = 7 (-u a + v a ). Here, v a is the projected relative 
velocity (in the u a frame), and 7 is the associated Lorentz factor. 

For a given photon momentum p a = E(u a + e a ), the energy E and direction of 
propagation e a in the u a frame are given by 

E = 1 E(l + e a v a ) (21) 
e a = [ 7 (1 + e^)]- 1 (e a + u a ) - 7 (« a + v a ). (22) 

The screen projection tensor H ab in the u a frame is related to the equivalent tensor 
in the u a frame by 

n ab = n ab - ^ P{a n b)c v c + (l^j PaPb n c ^v cl v C2 . (23) 

The polarization tensor in the u a frame, P a b(E, e), is related to P a b(E, e) by 

E- 3 P ab (E, e) = E~ 3 n c a ^P ClC2 (E, e). (24) 

From this transformation law we recover the frame-invariance of I(E, e) / E 3 , and the 
degrees of circular and linear polarization. 

It is straightforward to verify that under changes of frame the action of the 
Liouville operator transforms according to 

C[E- 3 P ab (E,e)\ = H c a mt 2 C[E- 3 P clC2 (E,e)}, (25) 

so that the collision tensor in the Boltzmann equation must transform as 

K ab (E,e)=H<?Ht 2 K clC2 (E,e). (26) 

This transformation law provides a useful means of generating the scattering tensor 
for a general choice of frame given its form in some particular frame. (It is often the 
case that the scattering process picks out a frame in which K ab (E, e) is particularly 
simple.) The transformation laws for the radiation multipoles Ia { (E), Va ; (E), Sa l (E) 
and Eai (E) are non-local in energy and show multipole coupling due to relativistic 
beaming. Since the detailed form of the multipole transformation laws is not needed 
for the calculation of linearised CMB anisotropies and polarization, we shall not give 
them here. 
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2.2.2 K ab (E,e) for Thomson Scattering 

Over the epoch of interest for the generation of CMB anisotropies and polarization, 
the dominant scattering mechanism for CMB photons is Compton scattering off free 
electrons. Well after electron/positron annihilation it is a good approximation to 
ignore electron recoil and the small thermal effects due to the electron distribution 
function. In this limit it is convenient to work initially in the rest frame of the 
tightly-coupled electron/baryon plasma, the four-velocity of which we denote by u a . 

The scattering tensor K ab (E, e) can be deduced from the standard result for the 
Thomson scattering kernel expressed in terms of the Stokes parameters (see, for ex- 
ample, |T7|]). Assuming an unpolarized distribution of electrons, we find the following 
expression for the scattering tensor in the u a frame in 1+3 covariant irreducible form: 



-I(E,e) + I{E) 



E 2 K ab (E,e) = n e a T <^-^H a b 

+ l/ CiC2 (E)e Cl g C2 + h clC2 (E)~e c ^ 



+ 



-^(^g) + ^[/ Q6 ^)] TT + ^[^(^)] TT 



2 ^abc& 



-V{E ) ~e) + \v d {Eu 



~d 



(27) 



where n e is the electron number density in the u a frame, and ot is the Thomson cross 
section. Equation ( p7|) is an exact expression for the scattering tensor in the rest 
frame of the electrons. The first set of terms in square brackets source the intensity 
evolution, the second set source the evolution of the linear polarization, while the 
final set source the circular polarization. Note that the source terms for the intensity 
have no monopole (Z = 0) moment in the u a frame, since there is no energy transfer 
in the Thomson limit in the electron rest frame. Note also how the quadrupole 
(I = 2) moment of the electric polarization sources the intensity, and the quadrupole 
moment of the intensity sources the electric polarization. (It is for this reason that 
polarization should be included in accurate CMB codes even if only the temperature 



information is of interest.) Equation fl27|) is equivalent to the scattering terms derived 
recently within the context of the "total angular momentum" advocated by Hu and 
collaborators 



19| , equation ( f27|) is free 



fig , pT9| . However, unlike the results in | 
from any (spatial) harmonic decomposition of the radiation variables. 

Transforming the scattering term to an arbitrary frame, and integrating over 
energies, we find the exact evolution of the intensity: 

J dE E 2 C[E- 3 I(E,e)] = -n e a T l(l + e a v a ) J dEI(E,e) 



+ — n e or T [ 7 (l + e\)]" 3 (/ + ( ClC2 e cl e c 



(2? 
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the linear polarization: 



J dEE 2 C[E- 3 V ab (E,e)) = -n e cr T y(l+e c v c ) J dEV ab (E,e) 



+ ^n e( T T [ 7 (l + e% c )}- 3 mn c b X ClC2 ] TT } (29) 

and the circular polarization: 

J dEE 2 C[E- 3 V(E,e)} = -n e a T 7(l + e a v a ) J dEV(E,e) 

- -jU e a T [ 7 (l + e a v a )]- 3 V b e b . (30) 

Here, tildes refer to quantities in the rest frame of the electrons, and 

3 ~ 9 ~ 

Cab = ~^ab + 2^ ab ' 

Expressions for /, I ab , and £ ab in terms of quantities in the u a frame can be derived 
from the multipole transformation laws if required. 

Equation (|30|) demonstrates that circular polarization is not sourced by Thomson 
scattering off unpolarized electrons, so that in the absence of polarising agents, such 
as primordial magnetic fields, we should not expect the primordial component of the 
microwave sky to have any circular polarization. In light of these remarks, we shall 
not consider the V Stokes parameter any further here. (Some effects of primordial 
magnetic fields on the CMB are discussed in [EO, pi]].) 



3 Linearisation Around a FRW Model 

The observed isotropy of the CMB (the anisotropies are at the level of only one 
part in 10 5 combines with the Copernican assumption to imply that on large 
scales the universe is well approximated by perturbing a FRW model [0, In the 
1+3 covariant and gauge- invariant approach to studying perturbations in cosmol- 



ogy |23|, £24], ^], the perturbed model is investigated in terms of covariant variables 
which are defined in the real universe (as opposed to the background model), and 
which vanish identically in the background model. With a physical choice for the 
fundamental velocity u a , the approach is gauge-invariant to all orders. Given a 
smallness parameter e for the almost-FRW model, the covariant variables that de- 
scribe inhomogeneity and anisotropy are at most 0(e). Important examples of 0(e) 
variables in a generic almost-FRW spacetime include the kinematic variables derived 



8 



from u a : 

shear o ab = D( a u b) 

vorticity u ab = D [a u b] 

acceleration A a = u a 

expansion inhomogeneity Z a = SD a Q ; (32) 

the matter variables derived from the matter stress-energy tensor T ab : 

momentum density q a = T( a ) b u b 

anisotropic stress n ab = T/ ab \ 

density inhomogeneity X a = SD a p/p; (33) 

the Weyl tensor C a bcd and hence its electric and magnetic parts: 

electric E ab = C acbd u c u d 

magnetic H ab = -\e acd C be cd u e . (34) 

Here, = V a w a is the fractional volume expansion rate, p = T ab u a u h is the mat- 
ter energy density, and S is a covariantly defined scale factor, S/S = 0/3 with 
D a S = 0(e). The linearised perturbation equations for the 0(e) variables are derived 
from the exact 1+3 covariant hydrodynamic and gravito-electromagnetic equations 
(see, for example, fli~2"|). The equation set closes when supplemented by constitutive 
relations for the matter variables. 



3.1 Linearised Multipole Equations 

The equations of radiative transfer provide the constitutive relations for the radia- 
tion sector of the cosmological model. It turns out to be most convenient to work 
directly with the evolution equations for the radiation angular multipoles. The ex- 
act multipole equations can be derived from equation ([18]). The equations for the 
intensity multipoles Ia 1 (E) are derived in 0, |10|, ; those for the polarization mul- 



tipoles are given in ||26|| . Here, we shall only make use of the linearised multipole 
equations. Since Thomson scattering does not produce any spectral distortion at 
linear order (this follows from equation (0)), it is sufficient to consider only the 
energy- integrated multipoles Ia { , £ai and 6i r Noting that Sjy x and Ba { are 0(e) 
for all / > 2 in an almost-FRW universe, while Ia { are 0(e) for I > 0, we find the 
following linearised evolution equations for the electric polarization^]: 

e * + 3 0Sa ' + ([Tip D £bAl " (27+T) D<(A - 1> 

curli3 A , = -n e a T \£ Al - ^zSnCab J , (35) 



2 We use the convention that Ea, and Ba, vanish for I < 2. 



9 



and for the magnetic polarization: 



! + 3°^ + (/Tip Mi " (27Ti) D<A - 1> 

2 

+ jj - jr curl £4, = -n e a T B Ar (36) 

Here, the curl of a rank-/ tensor is defined by 

cutISa, = e ClC2 ( aiJ D Cl S , j4i _ l) C2 , (37) 

which is PSTF for S Al PSTF. Equations (|35| ) and (|36| ) reveal that £a ; and B A{ 
are coupled through curl terms (as with Maxwell's equations), and that only the 
electric quadrupole £ a b has inhomogeneous scattering terms in its evolution equation. 
These observations are crucial for the prospect of detecting the imprint of primordial 
gravitational waves on the CMB (see Section |3.2| ) . 

The linearised intensity multipole equations are summarised by 

lAi + + D b I bAl - (2/ + jT -PfrA-i) + ^ IA aAl - J^I<7 ai aJl2 

= ~n e a T ^I Al - ^Iv ai 5 n - ^I( ai a 2 $i2^ , I > 1, (38) 

where v a is the projected relative velocity of the electron rest fame. The intensity 
monopole (radiation energy density) evolves according to 

i+iei + D a I a = 0, (39) 

which implies the following evolution equation for the comoving projected gradient 
of the energy density: 

(SDjy + ^QSDJ + SD a D b I b + pZ a - ^QSIA a = 0. (40) 



The linearised multipole equations (pq) , (p6|), (|38|) and fl40|) provide a complete 
description of linearised CMB physics in an almost-FRW universe. From these equa- 
tions one can infer the model-independent behaviour of the anisotropy and polariza- 
tion. Prior to recombination, the photon mean free time is short compared to the 
characteristic expansion time of the universe. On scales much larger than the pho- 
ton mean free path, the photons are tightly-coupled to the electron/baryon plasma, 
and the combined system behaves like a fluid. On super-horizon scales D a I grows 
due to gravitational instability, while on sub-horizon scales power oscillates between 
D a I and I a (acoustic oscillations). The fluid approximation breaks down inside the 



diffusion scale |p7| , where higher-/ intensity multipoles are excited at the expense 
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of damping the acoustic oscillations. During recombination the mean free path of 
the photons rises to effectively infinity, with the result that power is transported 
to higher-/ intensity multipoles through the free-streaming projection of D a I, the 
electron peculiar velocity v a , and the combination of intensity and electric polar- 
ization quadrupoles ( a b from the last scattering surface. Power is also transferred 
to the intensity multipoles through variations of redshift since last scattering along 
different lines of sight. The electric polarization quadrupole E a b is generated through 
recombination from the growing intensity quadrupole anisotropy. This mechanism 
ceases at last scattering, after which the power in the electric quadrupole is redis- 
tributed amongst the higher-/ electric and magnetic multipoles via free-streaming. 
Any reionisation tends to damp the primary intensity anisotropies while boosting 
the polarization signal, as well as generating new secondary intensity anisotropies 
[28|j for a recent review). 



sec 



3.2 Scalar Perturbations 

Up to this point our discussion has been concerned with arbitrary (small) perturba- 
tions around a FRW model. In particular, unlike most other perturbative approaches 
(for example, fl29fl), we have not had to split the perturbations into scalar, vector or 
tensor modes (see, for example, fl30| ) to derive the basic, gauge-invariant perturba- 
tion equations. However, for detailed calculation it is convenient to make a non-local 
mode decomposition of the 0(e) covariant variables to decouple the temporal and 
spatial dependencies of the 1+3 equations. The decomposition can be performed in 
linear theory in a fully covariant manner, which leads to a physically transparent, 
gauge-invariant perturbation theory. Here we shall only give a brief review of the 



procedure for scalar modes (those which describe density clumping) 0, [24j. Vector 
and tensor perturbations are treated in a similar manner. 

By definition, in scalar modes the 0(e) variables are constructed from 0(e) 
scalar potentials. It follows from the constraint equations of 1+3 covariant hy- 
drodynamics [T2j that the vorticity vanishes at 0(e) in a scalar mode. Further- 
more, since all curls vanish by construction, the constraint equations of 1+3 gravito- 



electromagnetism |12], |31| constrain the magnetic part of the Weyl tensor to vanish 
at linear order. The fact that curl£/i ; and curlf?^ vanish at linear order removes 
the coupling between the electric and magnetic polarization multipoles in their evo- 
lution equations (|35|) and ([36]). It follows that the hierarchy for the Ba { has no 
inhomogeneous terms present, so that scalar modes do not generate magnetic polar- 



ization 1 15, 16 1 . Any detection of primordial magnetic polarization would provide 



evidence of vector or tensor modes. This proves to be a much more promising ap- 



proach than relying on their cosmic variance limited temperature signature |32 



The scalar potentials from which the 0(e) covariant variables are derived may be 
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expanded in the complete set of scalar eigenfunctions Q" of the co moving Laplacian: 



D a D a Q™ = —Q 



S 2 



(41) 



with (j)^- 1 = 0(e). Rank-/ PSTF tensors, which are covariantly constant along the 
integral curves of u a , can be constructed from the Q^: 



Q%=(fy l D {ai ...D ai) QM 



(42) 



The Q Ai can be used to expand the radiation multipoles, so that, symbolically, 

e Al = iJ2£Ml (43) 

k k 

where the Ijp and are 0(e) scalars with 0(e 2 ) projected gradients. It proves 
convenient to use the notation lj® for the scalar coefficients in the expansion of 
D a I: D a I = I^2 k (k/ S)I k Qa , so that the mode-expanded intensity multipole 
equations take the form 



k 



f I (i-i) (l + l) 
S{(21 + 1) k (21 + 1) 



+ 3S Ak5n ~ 15 S^ 12 



rP+i) 



-n e a T [ ljp - ^v k 5 n - ^C k Si2 



2 

15 



(44) 



Here Aj., cr k , v k and ( k are the dimensionless coefficients in the mode expansions of 
A a , & a b, v a and Cafe, and 6K/S 2 is the zero-order 3-curvature scalar. The evolution 
equation for ij® follows from equation (pUf): 



k r(1 ) , 4 k 
35" 



4 0) + 4 --2, - -B.4, 0. 



4 
-( 

3 



(45) 



where i^, are the dimensionless coefficients in the expansion of Z a . For the electric 
polarization we find 



6(1) _ 



k 



I 



s{'(2l + l 



E 



_ (Z + 3)(i-l) 
(2/ + 1)(Z + 1) 



1-1(1 + 2)^ 



(i+i) | 



2 
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-n e a T ( 4° - ^( k 6 l2 ) . (46) 



At this point we have essentially recovered the well-known multipole equations 



for scalar perturbations in an almost-FRW model |19|, |33|. Arguably, the equations 
given here are somewhat more transparent than their gauge-dependent counterparts 
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10 100 1000 



Figure 1: Intensity Cf 1 and polarization Cf £ power spectra from scalar modes in 
a scale-invariant, adiabatic, open cold dark matter (CDM) model. The baryon 
fraction is = 0.05, the CDM fraction is Q c = 0.4, the Hubble constant is 
H = 50km s -1 Mpc -1 , and there is no cosmological constant nor reionisation. 

in alternative approaches. For example, the intensity monopole describes the 
spatial inhomogeneity of the radiation energy density, which is of clear physical sig- 
nificance, rather than the gauge- dependent difference between the density in the 
real universe and the background model. The mode-expanded multipole equations 
can be integrated directly to determine the anisotropy and polarization for some 
given initial conditions. To close the equations one must also integrate the mode- 
expanded hydrodynamical equations for the kinematic variables that source the in- 
tensity equations through redshift effects || . Since much of the complexity of the 
coupled multipole equations arises because they describe evolution along the time- 
like integral curves of u a , rather than along the lightlike trajectory of the radiation, 
it is more efficient to employ integral solutions to the multipole equations rather 
than integrate a large number of them directly [01 1. The integral solutions consist 



of simple source functions integrated against known geometric projection functions. 
The integral solutions for scalar and tensor modes are given in ||2?J; see also [F9 . 



An example of the intensity and polarization power spectra from scalar pertur- 
bations in an open model is given in Figure |l[ The code to produce the figure 
was based on CMBFAST p4]| ; the necessary modifications to solve the covariant 
equations were implemented by Antony Lewis. 
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4 Conclusion 



The 1+3 covariant formulation of radiative transfer provides a convenient, exact 
framework within which to study CMB physics. Although we have only considered 
the linearised form of the transfer equations in our discussion of the anisotropy and 
polarization in an almost-FRW universe, the formalism is of wider applicability. For 
example, the approach appears well suited to a systematic study of second-order 
geometric and scattering effects J7j. In the application to linearised CMB physics 
in an almost-FRW universe, the 1+3 covariant perturbative approach shares many 
of the advantages of the "total angular momentum" method |18|, [L9| . These include 
physical transparency, the explicitly quadrupole nature of the radiative source terms 
in the Boltzmann equation that arise from the angular and polarization dependence 
of the Thomson cross section, and the direct manner in which the angular multipoles 
describe the angular scale of anisotropy on the sky. However, the 1+3 covariant 
method goes further by ensuring manifest gauge-invariance, by deferring any non- 
local splitting of the perturbations until a late stage of the calculation, and by 
providing a straightforward linearisation procedure around a variety of background 
models. 
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